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STATISTICS OF VELOCITY FROM SPECTRAL DATA: MODIFIED VELOCITY CENTROIDS 
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ABSTRACT 

We address the problem of studying interstellar turbulence using spectral line data. We find a criterion 
when the velocity centroids may provide trustworthy velocity statistics. To enhance the scope of 
centroids applications, we construct a measure that we term "modified velocity centroids" (MVCs) 
and derive an analytical solution that relates the 2D spectra of the modified centroids with the 
underlying 3D velocity spectrum. We test our results using synthetic maps constructed with data 
obtained through simulations of compressible magnetohydrodynamical (MHD) turbulence. We show 
that the modified velocity centroids (MVCs) are complementary to the the Velocity Channel Analysis 
(VCA) technique. Employed together, they make determining of the velocity spectral index more 
reliable and for wider variety of astrophysical situations. 

Subject headings: turbulence - ISM: general, structure - MHD - radio lines: ISM. 
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1. INTRODUCTION 

The interstellar medium (ISM) is turbulent, as 
known from o bservations of non-thermal cloud velocities 
l)Larsonll992|) . and also expected on theoretical grounds, 
because of the very large Reynolds numbers (defined as 
the ratio of inertial to viscous forces). Understanding 
turbulence is crucial for understanding of many physi- 
cal processes, including energy injection into the ISM, 
non-photoelectric heating of the ISM, star formation, 
propaga tion of cosmic rays, and heat transport, (see re- 
view b y iVazauez-Semadeni et alJ feOOOlCho fc Lazarianl 
and references therein. 



Obtaining properties of ISM turbulence fro m observa- 
tions is a long-standing problem (see review bv lLazarianl 
1999, and references therein). The statistics of the ran- 
dom velocity field is essential for describing the turbu- 
lence. To obtain the ve locity statistics, so-called velocity 
centroids (Munch 1958) have been frequently attempted. 
However, the separation of the velocity and density con- 
tributions to the velocity centroids has always been a 
problem. Therefore the relation between the statistics 
of velocity and velocity centroids is frequently claimed 
to be trustwor thy only when the density fluctu ations are 
negligible fsee lVestuto. Ostriker fc Stonell2003j) . 

To remedy this situation a statistical technique 
termed Velocity C hannel Analysis, or in short VCA, 
was developed in iLazarian fc Pogosvanl (HXXJ. It 
has been successfully used to obtain the spectra 
of turbulence in HI in Small Magellanic Cloud 
and the Milky Way (see [S tanimir ovic fc Lazarianl 
l2001lDeshpande. Dwarakanath. fc Gossl 12000^ How- 
ever, the VCA requires t urbulence to be supe rsonic and 
to follow power laws (see lEsauivel et alll20f)8l ). 

This paper will identify a technique that can be reliably 
used with turbulence even if the turbulence is subsonic 
and/or does not obey a power-law. The latter case can 
occur when, for instance, self-gravity modifies the turbu- 
lence on small scales. In pursuing this goal we introduce 
modified velocity centroids (MVCs) and derive an ana- 

Astronomy Department, University of Wisconsin-Madison, 475 
N. Charter St., Madison, WI 53706, USA. 
lazarian@astro.wisc.edu; esquivel@astro.wisc.edu 



lytical relation between the statistics of MVCs and the 
underlying statistics of the 3D turbulent velocity. 

2. MODEL: TURBULENCE IN MAGNETIZED GAS 

There has been substantial progress in understanding 
of compressible magne tohydrodynamical (MH D) turbu- 
lence (see a review bv lCho fc Lazari an 2003a) that will 
guide us in adopting the model for this paper. 

Num erical simulations in iCho fc Lazarianl (2002, 
2003b) revealed anisotropic spectra of Alfven and slow 
waves, as well as an isotropic spectrum of the fast waves. 
The anisotropies are elongated along the local direction 
of magnetic field. However, since this local direction 
changes from one place to another, the anisotropy in 
the system of reference related to the mean magnetic 
field is rather modest, i.e. of the order of SB/B, where 
SB is the dispersion of the random m agnetic field. For 
typical ISM SB ~ B <|SDoelstralll972f) and the expected 
anisotropy is modest. Therefore it is possible to charac- 
te rize the turbulence us ing isotropic statistics (see testing 
in lEsauivel et alj f2003 and in §4). 

Only z-components of the velocity are available 
through Doppler shift measurements. In what follows 
we shall assume that the velocity and density of the gas 
can be presented as sums of a mean value and a fluctu- 
ating part: V z = v + v, p to t = Po + P, where v = (V z ), 
Po = (ptot)- The fluctuating components v and p satisfy 
the condition (v) = 0, (p) = 0. Henceforth we use (...) to 
denote ensemble averaging. The correlation function of 
the z-component of velocities at two points Xi and x 2 , 
separated a distance r (with r = xi — x 2 ), denoted by 
subscripts 1 and 2 (V Zy \ = vq + v±, V z ,2 — vq + V2) is 
(V1V2) = B zz (r). And it is related to the component 
of the 3D spectrum th rough a Fourier transform (see 
IMonm fc YatloTr1IT97^ : 



e lkr F zz (k)dk, 



(1) 



B zz (r) 

where 

F zz (k) = (F LL (k) 

is the projection of the spectral tensor, expressed through 
the transversal Fnn and longitudinal Fll (normal and 



F NN (k))k 2 Jk 2 + F NN (k) (2) 
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parallel to r, respectively) 3D spectra of the velocity field, 
which arc functions of the wavenumber amplitude k. 

Similarly, the two point density statistics can be given 
by the structure function ((pi — P2) 2 ) = d(r). 

3. MODIFIED VELOCITY CENTROIDS (MVCS) 

In what follows we consider the emissivity of our me- 
dia proportional to the first power of density 2 and no 
absorption present. In this case, the intensity, /(X) = 
/ p s dV , is proportional to the column density, where p s 
is the density of emit ters in the Position-Position Veloc- 
ity (PPV) space fsee lLazarian fc Pogosva n 2000). That 
is, /(X) = J p s dV — J aptotdz, where, ptot is the den- 
sity of atoms, a is a proportionality constant, dz denotes 
integration along the line of sight, and X is a two dimen- 
sional vector in the plane of the sky. 

The structure function of MVCs consists of two terms. 
It will be clear from the discussion below that it may be 
safe to use ordinary centroids when the second term is 
much smaller than the first term. In particular, 



M(R) = ( [5(X0 - S(X 2 )] 2 - (V 2 ) [J(X0 - J(X a ) 



( 3 ) 

where R = Xi — X2, S(X.) is a "unnormalized cen- 
troid" of the z-component of velocity V z , and (V 2 ) = 
v a + v Lrb + v thermai- To simplify our notations we 
shall omit the turb and z-subscripts (also for eq.©). 
We do not discuss the effect of vthermai separately be- 
cause for the technique presented here it acts in the same 
way as regular motion, and one can formally introduce 
v' 2 = v 2 + v 2 thermal . S(K) = J x Vp s dV, that is, at a 
given position X, integrating along the line of sight 3 . 

From an observational standpoint, the velocity disper- 
sion (V 2 ) can be obtained using the second moment of 
the spectral lines: 



(V 2 ) 



{J^V 2 p s dV)/{ 



PsdV). 



(4) 



To express the 2D statistics of the modified centroids 
M(R) through the underlying 3D velocity and density 
statistics we shall have to make several elementary trans- 
formations. First of all, the difference of integrals that 
enter the expression for the modified centroids should be 
presented as a double integral. Then using an elemen- 
tary identity (a-b)(c-d) = l/2\(a - d) 2 + (b - c) 2 - (a - 
c) 2 — (6 — d) } it is possible to find fsee lLazaria 5 11991 

((S(Xi)-S(X a )) 2 ) =a 2 J J dz 1 dz 2 [D(r)-D(\z 1 -z 2 \)}, 

(5) 

where D = ((Vip to t,i - VzPtnt.i) 2 )- Using th e Mil- 
lionshikov hypothesis (see lMonin & Yagloml ll975'l to re- 
late the fourth moments of the fields with the second 
moments, namely that {hxh^h^hi) « (hxhz) (/13/ki) + 

2 This is the case for the H I 21cm line. For recombination radi- 
ation (such as Ho for instance) the emissivity is proportional to the 
square of the density, and modifications of the present technique 
are necessary. 

3 In terms of practical data handling, to minimize the effect of 
absorption at the centers of the lines, subtracting of mean values 
of velocity vq in the expression for S(X) and (V 2 ) may be advan- 
tageous. In terms of S(X) the use of <5S(X) = 5(X) — (5(X)} may 
be preferable. 



(hxh3)(h2hi) + {h\hi)(h2h^) we get sufficient simplifi- 
cations of the expression: 

D « <y 2 )d(r) + (p 2 ot }((^i-^) 2 )-i((yi-^ 2 ) 2 )d(r)+ C (r), 

(6) 

where c(r) denote the cross terms arising from correla- 
tions between the fluctuations of the velocity and density 

c(r) = 2{vip 2 ) 2 - 4w (wipip 2 ) - 4po(piWiw 2 ). (7) 

The correlations between velocit y and density fluctu- 
ations have been st udied earlier ijEsauivel et alJ 120031: 
lLazarian et al.ll2001() and found not to have a strong im- 
pact on the VCA. In parti cular, the first term in eq. 
Q) was directly obtained in lEsauivel et al.1 {2003), and 
numerical integration of the corresponding fluctuations 
provides fluctuations in the centroid value that does not 
exceed 8% of the MVCs values at small \R\. In this paper 
we will neglect the cross terms and test this assumption 
by analyzing synthetic maps of compressible MHD simu- 
lations. A more detailed analysis, with more sets of data, 
and the cross terms, will be included in a forthcoming 
paper. 

Eq. 10 depends on both density and velocity statistics. 
The dominant contribution arising from density fluctua- 
tions is the first term of eq. ©. Integration of this term 
yields the second term in eq. © , as result the expression 
for the MVCs is 



M(R) » a 2 (V 2 }(p 2 tot ) / / dz,dz 2 [b(r) ~ b(\ Zl 



z 2\ 



where 



6(r) = 



((V^V 2 ) 2 ) 



1 _ I d{r) 



(8) 



(9) 



(v 2 ) V 2 ( P | ot ) 

If {l/2)(d(r)/(p 2 ot )) < 1 then eq.® recovers the statis- 
tics of the velocity field. This is certainly the case if the 
turbulence is subsonic, where the amplitude of the den- 
sity fluctuations is small, or when the density spectrum 
is steep (i.e. most of the power residing at large scales) 
given that we are interested in the fluctuations at small 
R. 4 However, if the amplitude fluctuations of density 
at small scales is large (i.e. shallow spectrum or highly 
supersonic turbulence) , the contribution could be impor- 
tant. If we omit the density term in eq.® and take a 
2D Fourier transform to resulting distribution of MVCs 
we have: 



P 2 (K) = J d 2 R e- iK R M(R), 



(10) 



where we have used again the notation of using capital 
letters to distinguish 2D vectors from the 3D counter- 
parts. Substituting eqs. JSJ, (JTJ, and into eq. JTHJ 
one gets 

P 2 {K) ~ C 1 F NN (K,0) + C 2 5(K) (11) 

where C\ and C 2 are constants (C\ — 4na 2 (p 2 ot )L, 
where L is the extent of the emitting region) and the 
second term provides only a 5 function contribution 
at K = 0. Thus eq. Ijlljl provides a relation 5 be- 
tween the spectrum of the modified centroids and the 

4 By steep spectrum we m ean P30 oc k n with n < —3. If n > —3 
the spectrum is shallow (see Lazarian & Poeosvan 2000). 

5 Our result is similar to that in Kitamura 1 1993) who used 
normalized centroi ds given bv eq. j l2| but implicitly assumed that 
/' 



: const (see also Vcstuto et al. 2003). 
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transversal spectrum of the 3D velocity field. We should 
note that the relation in eq. I jlljl does not depend 
on a particular form (power-law is a common assump- 
tion) of the underlying spectra. If the turbulent yeloc- 

~~ H[l996; 



ity field is mostly solenoidal (see |M atth; 
iPorter. Woodward fc Pouauetlll998t) . the isotropic spec- 
trum E(k) that is used for describing hydrodynamic tur- 
bulence is uniquely defined through F^jq(k), namely, 
E(k)^Ank 2 F NN {k). 

4. NUMERICAL TESTING 

In what follows we will in clude the normalized cen- 
troids in their usual form fsee lMiesch fc Ballvlll994() 



C(X) 



V Ps dV/ / Ps dV. 



(12) 



Intuitively it is clear that the power spectrum of the cen- 
troids given by eq. I|12|) should provide a better fit to the 
power spectrum of the velocity (the spectral index, not 
the amplitude) than that of <S(X). Indeed, the contribu- 
tions of density fluctuations are mitigated by the division 
over the column density in eq. (|12|l . 

The testing was done using the files obtained through 
simulations of compressible MHD turbulence. The de- 
scription o f the code and the num erical sim ulations can 
be fou nd in lCho fc Lazarianl l|2002|) and in lEsauivel et al.l 
(12003ft . We u s e the same data cube that we used in 
lEsauivel et al.l l)2003|) . Namely, the dimension of the 
cube is 216 3 , and the Mach number 2.5. The result 
of the MHD simulations are velocity, density and mag- 
netic field data. We use density and velocity to produce 
spectral line data cu bes, the procedure is described in 
lEsauivel etaDpOOl . 

We calculated 5(X), C(X), and /(X) from synthetic 
spectral line data cubes. The spectra -P2(K) were ob- 
tained by applying a fast Fourier transform to the 2D dis- 
tribution of the centroids and column density values. The 
power spectrum of the MVCs was calculated subtracting 
(V 2 ) times the column density power spectrum, from the 
power spectrum of 5(X) (see eq. |3J). For the original 
data cubes (with a steep density spectrum), the differ- 
ences among all three types of centroids were marginal 
for the dynamical range studied, although unnormalized 
centroids systematically show more deviations from the 
velocity power spectrum. To enhance the effect of den- 
sity, we made the density shallow by changing the Fourier 
components of the density amplitudes, but, in order to 
preserve the density-v elocity correlations, we kept the 
phase information fsee lEsauivel et aLll2003(l . The results 
are shown in Fig 1. It is clear that the three measures 
provide a good fit at large scales, but at small scales 
normalized and unnormalized centroids tend to be den- 
sity dominated. An additional advantage of using MVC 
over the conventional normalized centroids is that MVC 
provides the correct amplitude of the spectrum, while 
normalized centroids at best give the correct index (log- 
arithmic slope). 

It is also clear that for a shallow density spectrum, the 
improvement from unnormalized centroids to normalized 
centroids is marginal. 

Our numerical study shows that MVCs are advanta- 
geous when the underlying spectrum of density is dom- 
inated by fluctuations at small scales. And, the main 
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Fig. 1 . — The comparison of the 2D spectrum of the three types of 
centroids with the underlying 3D spectrum of velocity fluctuations. 
The solid line corresponds to the power spectrum of velocity for 
compressible MHD simulations with Mach number ~ 2.5. The 
dotted line corresponds to the spectrum of density modified by 
changing the amplitudes the Fourier components of the original 
spectrum. The modified density spectrum E(k)/k 2 ~ fc~ 2,5 for 
the inertial range. The vertical position of the normalized centroids 
(triangles) is shifted vertically for visual purposes. 



advantage of MVCs is that they allow an explicit statis- 
tical description of the procedures involved, which makes 
the analysis more reliable. 

5. DISCUSSION 

Our study revealed that the MVCs can be successfully 
used to obtain the velocity statistics from synthetic ob- 
servational data. Another technique, namely, the VCA 
is complementary to the MVCs. VCA is more robust to 
de nsity-velocity correlati o ns. Fo r instance, it was shown 
in lLazarian fc Pogosvanl l)2f)f) f)) that VCA can provide 
correct statistics of velocity fluctuations even in the case 
when the velocity and density are correlated at their 
maximum given by the Cauchy-Swartz inequality (see 
iMathews fc Walkerll97Tllh On the contrary, if we assume 
such a maximum correlation for the MVCs we will not 
get the correct scaling. Comparing results obtained with 
the MVCs and the VCA one can get a better handle on 
what are the velocity statistics. 

It is advantageous to combine VCA and MVCs tech- 
niques. We can apply VCA to the largest scales, where 
turbulence is supersonic, and study MVCs at all sca- 
less, including those where the turbulence becomes sub- 
sonic. The correspondence of the spectral slope obtained 
with the two different techniques would substantially in- 
crease the reliability of the result. There is a more sub- 
tle, but important point. MVCs provide the spectrum of 
solenoidal motions given by Fnn> while VCA is sensitive 
to both potential and solenoidal motions. This may be 
used to get both potential and solenoidal motions in or- 
der to estimate the role of compressibility. In addition, 
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because the thermal velocity acts in the VCA analysis 
the same way as the thickness of the channel maps, mea- 
surements of the spectrum with the MVCs can provide a 
means of estimating the thermal velocity. It is essential 
to combine the advantages of the techniques. 

Our analysis of numerical data in § 4 shows that the 
worries about the velocity centroids ( see iStenholmll 1 990(l 
although justified, may be somewhat exaggerated. 

If we assume that the criterion for the use of the cen- 
troids is satisfied (that the second term in eq. 10 is 
small compared to the first one) for earlier data sets, the 
analysis of observational data in ijMiesch fc Ballvlll994() 
provides a range of power-law indexes. Their results ob- 
tained with structure functions if translated into spec- 
tra are consistent with E{k) = fc' 3 , where [3 = —1.86 
with a standard deviation of 0.3. The Kolmogorov in- 
dex —5/3 falls into the range of the measured values. 
L1228 exhibits exactly the Kolmogorov index —1.66 as 
the mean value, while other low mass star forming re- 
gions L1551 and HH83 exhibit indexes close to those of 
shocks, i.e. ~ —2. The giant molecular cloud regions 
sho w shallow indexes in the ran ge of —1.9 < (3 < —1.3 
CseelMiesch. Scalo fc Ballvlll999D . It worth noting that 
iMiesch fe^Beu4vll)1994|) obtained somewhat more shallow 
indexes that are closer to the Kolmogorov value using 
autocorrelation functions. Those may be closer to the 
truth as in the presence of absorption in the center of 
lines, minimizing the regular velocity used for individ- 
ual centroids might make the results more reliable (see 
also first footnote in §3). Repeating of the study using 
MVCs and taking care of the absorption effects may be 
advantageous. 

The emergence of the Kolmogorov index for the com- 
pressible magnetized gas may not be so surprising. 
iCho fc Lazarianl pOO^ for a range of Mach numbers re- 
ported the existence of distinct Alfvenic, slow, and fast 



wave cascades. The Alfven and slow waves result in 
Kolmogorov-type scaling. The application of the VCA 
to S mall Magellanic Cloud iStanimirovic fc Lazarianl 
|2001[) resulted in the velocity spectrum consistent with 
those predictions. However, other researchers (e.g. 
iVestuto et alJ l2003j) reported indexes that varied with 
the Mach number of the simulations. It is clear that 
more theoretical and observational work is necessary. 

6. SUMMARY 

1. We derived a criterion when the velocity centroids 
may reflect the actual underlying velocity statistics. We 
introduced MVCs that may recover velocity statistics 
from spectral line data in cases when the traditional cen- 
troid analysis fails. 

2. Our numerical tests show that both the MVCs and 
the normalized velocity centroids (see eq. C2D success- 
fully recover the velocity information from data cubes 
obtained via numerical simulations of the compressible 
MHD turbulence if the density spectrum is steep. 

3. MVCs and the VCA are complementary indepen- 
dent techniques. Obtaining the same power spectrum 
with both of them enhances confidence in the result. 
MVCs are most useful when turbulence is subsonic 
or/and does not follow the power law. Additional 
information on the solenoidal versus potential motions, 
temperatures of gas etc. may be obtained combining 
these two techniques. 
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